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I present the results of searches for B meson decays into two charmless mesons. I take 
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into account final states made up of two pseudo-scalar (PP) or one pseudo-scalar and 
one vector (PV). The measurements use the data samples collected and analysed by the 



B-factory experiments at the T(4S) resonance energy: BABAR, Belle and CLEO. 
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^ ■ 1 Introduction 

^| Charmless hadronic final states play an important role in the study of CP-violation. In the 
Standard Model, all CP-violating phenomena are a consequence of a single complex phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1. 

Charmless decays of B mesons may proceed by b — > u, b — > s, or b — > d transitions. The 
latter two mechanisms require flavor changing neutral currents which are not present at tree 
level in the Standard Model, and therefore must occur through suppressed processes such as the 
penguin mechanism. Such processes involve loops, which can get contributions from physics 
beyond the Standard Model. Even in the absence of such new physics, interference among 
competing amplitudes for a given decay mode can be exploited to measure CKM phases. 

These charmless B decays in which CKM favored amplitudes are suppressed or forbidden 
are rare decays and they typically have branching fractions (£>) of less than 10~ 4 . Even if the 
branching fraction for these modes is expected to be so low, the present data samples available 
to BABAR and Belle Collaborations together with the new analysis effort by CLEO allow for 
measurements or stringent limits on many such modes. 

In these years, the BABAR and Belle collaborations have reached solid results j2| on mea- 
surements of CP-violating asymmetries in B decays into final states containing charmonium, 
leading to strict constraints on the angle (3 of the CKM Unitarity Triangle. 

In the Standard Model picture, the angle a can be related to time-dependent CP-violating 
asymmetries in the analysis of the decay B° — > tt + tt~ as well as in the case of three body 
7T + 7r~7r° decays in which such measurement would exploit interference between the B° — > p ±r K T 
modes and the color-suppressed B° — > p°7i°. 
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Ratios of branching fractions for various tttt and Kir decay modes are in principle sensitive 
to the angle 7 j^j, even if the measured branching fractions of all the Kir modes show a clear 
path which seems to suggest the dominance of the penguin amplitude over the others As 
a matter of fact the Kn modes (together with the relative PV modes like K*ir, pK, ujK) are 
sensitive to 7 through the Cabibbo-suppressed amplitude term which is proportional to V* b V us 
but an interference between the latter and the Cabibbo-allowed amplitude term is needed in 
order to extract 7 1 . 

The extraction of a from measurements of the time-dependent asymmetry in B° — ► 7r + 7r~ 
is complicated by the interference of tree and penguin amplitudes with different weak phases. 
The B — > KK decays are characterized by similar penguin processes and, hence, can be used to 
isolate the tree and penguin contributions to B° — > 7i + tt~ . For example K°K is a pure b — > d 
penguin. The fact that no tree amplitude is possible, strongly reduces the expected branching 
fraction, but also allows to estimate the size of penguin contributions in B — > tttt, once the 
B(B — > K°K ) has been measured and SU(3) breaking effects have been taken into account. 

has the same penguin as K°K plus an annihilation contribution. The comparison of the 
two branching fractions can be used to understand which is the real contribution coming from 
annihilation terms. Moreover, the entire K !^ amplitude is equal to the penguin contribution 
to the 7T 7r° channel, assuming SU(3). So, in the scenario of a large branching fraction for 7r 7r°, 
the measurement of the B(K°K°) can help to disentangle penguin contribution from other 
subleading terms in 7r°7r° amplitude. These other terms can be deduced in a similar way from 
K + K~ decay which is a pure annihilation mode 

Also the decays of I?-mesons to pseudo-scalar and vector particles provide opportunities for 
investigating the phenomenon of CP violation. Of particular interest are the quasi-two-body 
B — ► pit decays to three-pion final states. These can be used to measure the angles a and 7. 

a can be determined from a full Dalitz plot analysis of the modes B° — > p ± 7r T and B° — > 
p°7T jH]. The data sets available are still too small for this type of analysis. For the time 
being though, the BABAR experiment also performed a simultaneous measurement of branching 
fractions and CP- violating asymmetries in the decays B° — > p^^ and B° — > p~K + (and their 
charge conjugate) using a time-dependent maximum likelihood analysis. Following a quasi-two- 
body approach [7j, the BABAR technique is restricting the analysis to the two regions of the 
7r =F 7r°/z. =t: Dalitz plot (h = it or K) that are dominated by either p + h~ or p~h + . From this kind 
of analysis, the time-integrated and flavor-integrated charge asymmetries A P J P and A p ^p that 
measure direct CP violation can be extracted together with the quantities S pn and C pn for the 
p7r mode. 

Since, like in the two-body case, the extraction of a from p ± n T is complicated by the 
interference of decay amplitudes with different weak and strong phases, another strategy is to 
perform an SU(2) analysis using all pn final states !9J. Assuming isospin symmetry, the angle a 
can be determined free of hadronic uncertainties from a pentagon relation formed in a complex 
plane by the five decay amplitudes B° — > p + 7r~, B° — > p~7r + , B° — * p°7r°, B + — > p + 7r° and 
B + — > p°7r + . These amplitudes can be determined from measurements of the corresponding 
decay rates and CP-asymmetries. 7 can be extracted from the interference of B + — > Xco 7r+ an d 
B + — > p°7r + in B + — > 7r + 7r~7r + decays 

1 For an example on this amplitude formalism see [S] 
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Also measurements of the charge asymmetries in the decay rates are reported according to 
the definition: 

•Acp — — 7= =\ ; ■ 

Here and throughout this paper charge conjugate modes are implied. We also make use of 
the notation h to represent a charged hadron that may be either a kaon or pion. 



2 Analyses 

The measurements presented in this paper are mainly based on the following data samples. 
The analyses of BABAR are based on (87.9 ± 1.0) x 10 6 BB pairs collected between 1999 and 
2002 with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. This 
corresponds to an integrated luminosity of approximately 81 fb -1 at the T(4S) resonance. The 
results from Belle are based on 78 fb" 1 data sample which corresponds to 85 x 10 6 BB pairs. The 
Belle detector ^T] is a large-solid- angle general purpose spectrometer at the KEKB asymmetric- 
energy e + e~ storage ring. The current CLEO data set is made up from the full CLEO II and 
CLEO III data samples totaling 15.3 fb -1 . CLEO[12j is a general purpose solenoidal magnet 
detector operated at the Cornell Electron Storage Ring (CESR), a symmetric-energy storage 
ring tuned to provide center of mass energies near the T(4S) resonance. 

Hadronic events are selected based on track multiplicity and event topology. Tracks are 
identified as pions or kaons using the various Particle Identification (PID) techniques described 
below. Candidate K® mesons are reconstructed from pairs of oppositely charged tracks that 
form a well-measured vertex and have an invariant mass within a defined window around 
the nominal K® mass ^H]- The candidate must have a displaced vertex and flight direction 
consistent with a K® originating from the interaction point. Candidate ir° mesons are formed 
from pairs of photons with an invariant mass within a defined window around the nominal 
7T° mass. The tt° candidates are then kinematically fitted with their mass constrained to the 
nominal ir° mass j!6j . 

Reconstructed B candidates are identified using two kinematic variables: the beam-energy 
substituted mass (or beam-energy constrained mass) titles = J (-^beam) 2 ~ (Pb M ) 2 aric ^ the en ~ 
ergy difference AE = Eg™ — E™ m , where E™ m is the beam energy, p§ M and Eg M are the 
momentum and the energy of the reconstructed B meson in the center-of-mass system (CM). 

In these charmless modes, the largest source of background comes from random combina- 
tions of tracks and neutrals produced in the e + e~ — > qq continuum (where q — u, d, s or c). 
These backgrounds are suppressed using the event topology. In the CM frame this background 
typically exhibits a two-jet structure that can produce two high momentum, nearly back-to-back 
particles, in contrast to the spherically symmetric nature of the low momentum T(4S) — > bb 
events. This topology difference is exploited by making use of event-shape quantities. Various 
techniques and variables are used. 

In the BABAR analyses, a preliminary requirement is applied on the angle 6 sp h ^3] between 
the sphericity axes, evaluated in the CM frame, of the B candidate and the remaining tracks 
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Figure 1: Examples of plots from the various two-body analyses and techniques. See text in 
Sect. E]for details. 



and photons in the event. The distribution of the absolute value of cos6 sph is strongly peaked 
near 1 for continuum events and is approximately uniform for BB events. 

Another similarly powerful variable is the flight direction of the B candidate given by 
cos 6b = p ■ z where p is the vector sum of the daughter momenta and z is the beam axis. 
Since the vector T(4S) is produced in e + e~ annihilation it has a polarization J z = ±1, and 
the subsequent flight direction of the pseudo-scalar B mesons is distributed as |Y 1 ±1 (6', 0)| 2 ~ 
sin 2 8 = 1 — cos 2 9. Background events are flat in this variable. 

Another quantity largely used in the charmless analyses is a Fisher discriminant T which 
consists of an optimized linear combination of variables that distinguish signal from back- 
ground [PTj . 

BABAR analyses use a two- variable Fisher discriminant with Y^iPi an d J2iPi \ cos^l 2 where 
Pi is the momentum and 9i is the angle with respect to the thrust axis of the B candidate, both 
in the CM frame, for all tracks and neutral clusters not used to reconstruct the B meson. 

Belle's technique consists of forming signal and background likelihood functions, Cs and 
Cbg-i from two variables. One is a Fisher discriminant determined from six modified Fox- 
Wolfram moments [13] and the other is cos8b- Requirements are imposed on the likelihood 
ratio LR = £ s /(£s + £-bg) f° r candidate events. 

CLEO builds its Fisher discriminant as a linear combination of fourteen variables including 
the direction of the thrust axis of the candidate with respect to the beam axis, cos8 t h r , and 
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mode 




CLEO 


BR (10~ 6 ) [UL 
BABAR 


@ 90% CL] 
Belle 


WA 


£? u ^ 


7T + 7r 


A c+1.4+0.5 
^•<J_1 2-0 4 


4.7 ±0.6 ±0.2 


4.4 ±0.6 ±0.3 


4.6 ±0.4 


B"^ 


K + ti- 


18.U_2 ]_ g 


17.9 ±0.9 ±0.7 


18.5 ± 1.0 ±0.7 


18.2 ±0.8 


B° -> 


K + K- 


< 0.8 


< 0.6 


< 0.7 


< 0.6 


B+ - 


K°7T + 


1 8 8+ 3 -' / + 2 - i 
J-O-O-3.3-1.8 


22.3 ±1.7 ±1.1 


22.0 ±1.9 ±1.1 


21.8 ±1.4 


fi+ - 


K + 7T° 


10 Q +2.4+1.2 
lZ.y_ 2 . 2-1.1 


12.8l^±1.0 


12.8 ±1.4li;3 


12.8 ±1.1 


£? u ^ 


K°7r° 


i o Q+4.U+1.7 
iz -°-3.3-1.4 


11.4 ±1.7 ±0.8 


12.6 ±2.4 ±1.4 


11.9 ± 1.5 


B + - 


K°K+ 


< 3.3 


< 2.5 


< 3.4 


< 2.5 


B ^ 


K°Tf 


< 3.3 


< 1.8 


< 3.2 


< 1.8 


B + - 


7T+7T 


/( £+1.8+0.6 
4 -°-l. 6-0.7 


5.5^ ±0.6 


5.3 ± 1.3 ±0.5 


5.3 ±0.8 


B ^ 


ttV 


< 4.4 


< 3.6 


< 4.4 


< 3.6 



Table 1: Branching Fraction results for PP charmless modes: tttt, Ktt and KK with all charge 
combinations JH] • Also the world averages are given. 



the nine conical bins of a "Virtual Calorimeter". They are the scalar sum of the momenta 
of all tracks and photons (excluding the B candidate daughters) flowing into nine concentric 
cones centered on the thrust axis of the B candidate, in the CM frame. Each cone subtends 
an angle of 10° and is folded to combine the forward and backward intervals. In addition, the 
momentum of the highest momentum electron, muon, kaon, and proton are used as inputs to 
the Fisher discriminant taking advantage of the high quality particle identification in CLEO 
III. 

In the case of a candidate mode involving one or more charged pions or kaons, such as 
B — > Kit or B — > tt-k , each charged track must be positively identified as K or rr. 

PID in BABAR analyses is accomplished with the Cherenkov angle measurement from a 
detector of internally reflected Cherenkov light (DIRC). The final fit to the data includes the 
normalized Cherenkov residuals (9 C — 9*) /(Jq c and [9 c — 9f\ /cro c , where 8 C is the measured 
Cherenkov angle of the charged primary daughter, <jg c is its error, and 9* (9^) is the expected 
value for a pion(kaon). The latter two quantities are measured separately for negatively and 
positively charged pions and kaons, from a sample of D° —>■ K~n + originating from D* + decays. 

PID in Belle experiment is based on the light yield in the aerogel Cherenkov counter (ACC) 
and dE/dx measurements on the central drift chamber. For each hypothesis (K or rr), the 
dE/dx and ACC probability density functions are combined to form likelihoods, Ck and C w . 
K or rr mesons are distinguished by requirements on the likelihood ratio £k/(£k + £n)- 

The K/tt identification in CLEO relies on the pattern of Cherenkov photon hits in the 
RICH detector fitted to both a kaon and pion hypothesis, each with its own likelihood Ck 
and C^. Calibrated dE/dx information from the drift chamber is used to compute a x 2 f° r 
kaon and pion hypotheses. The RICH and dE/dx results are combined to form an effective 
X 2 difference, A Kn = —2\ia£ K + 2\n C n + x\ ~~ xl- Kaons are identified by A Kn < 5k 
and pions by — A/^ < S^, with values of 5k and 5 n chosen to yield (90 ±3)% efficiency 
as determined in an independent study of tagged kaons and pions obtained from the decay 
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mode 




CLEO 


BR (10~ 6 ) [UL (6 
BABAR 


i 90% CL] 
Belle 


WA 


B ^ 




27.6+51 ± 4.2 


22.6 ± 1.8 ±2.2 


9D Q+ti.U+2.8 

Z<U. O a q q i 
— D.o — o. 1 


22.7 ±2.5 


B+ -> 


p U 7T + 


io.4^;2 ±2.1 


24 ± 8 ± 3 


8.01^ ±0.7 


9.5 ±2.0 


B+ -» 


P+7T 


< 43 






< 43 


B ^ 


p 7T° 


< 5.5 


< 10.6 


< 5.3 


< 5.3 


B ^ 


P+if- 


16t|±3 


7.3ti;l ± 1.2 


< 23 


7.7 ± 1.7 


B+ -» 


P °ir+ 


< 17 


< 29 


< 12 


< 12 


B+ -» 


p+if° 


< 48 






< 48 


B ^ 


p°K° 


< 39 




< 12 


< 12 



Table 2: Branching Fraction results from the ph modes with h being a n or a K with all charge 
combinations. 



D *+ ^ tt+^O pO ^ X-^). 

Finally PV modes allow for the use of the vector particle helicity angle 9h defined as the 
angle between the direction of one of h daughters in the h rest frame and the direction of h in 
the B rest frame. 

3 Results 

In general the charmless decay modes taken into account here have contributions from (a) 
signal, (b) continuum qq background and (c) cross-feed from other B modes. 

In order to extract the signal yields, BABAR and CLEO use unbinned extended maximum 
likelihood fit. The input variables to the fit are in general mESi AE, Fisher discriminant T . 
Some BABAR analyses take advantage of a Neural Network to fight the continuum background 
instead of the Fisher discriminant. In CLEO analyses, the flight direction of the B candidate 
cos#b is added to the fit inputs. Moreover most of the charmless decay analyses in BABAR 
include also the Cherenkov angle in the fit to distinguish when necessary among differences in 
the K/tt content of the final states. 

In Belle analyses, the signal yields are extracted by a binned maximum likelihood fit to the 
AE distribution in the ttles signal window (5.271 GeV/c 2 < m^s < 5.289 GeV/c 2 , the lower 
limit being 5.270 GeV/c 2 in case of a 7r° in the final state). The rriES distributions are fitted as 
a consistency check. 

Tables ^ |U and El show the Branching Fraction results of respectively the charmless two- 
body hadronic B decays [TH], the modes with a p in the final states [TH], the modes with a K*, 
a <p or an uj in the final states [2U|. The results shown are from all the three experiments and 
also the relative world averages are reported. 

In Figure Q plots are shown from the various two-body analyses and techniques. The 
top right left are two examples from BABAR analyses and represent the projections of the 
tues distributions from the unbinned maximum likelihood fit for events that satisfy optimized 
requirements on probability ratios for signal to background based on all variables except tties 
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mode 






BR (10~ 6 ) [UL (5 


) 90% CL] 








CLEO 


BABAR 


Belle 


WA 


B + - 




± 1.6 


15.5 ±3.4 ±1.8 


19.4+™1 


12.3 ±2.6 


B + - 




< 31 






< 31 


B°^ 


K* + TT~ 


16+^ ± 2 




< 30 


16 ±6 


B ^ 


K*°ir° 


< 3.6 




< 7 


< 3.6 


B+ - 


<pK + 


5.5±^±0.6 


10.0i^±0.5 


9.4 ± 1.1 ±0.7 


9.3 ±0.8 


B ^ 


(j)K° 


5.4+ 3 - 7 -2.7±0.7 


7.6±U ± 0.5 


9.0 ±2.2 ±0.7 


7.7 ± 1.1 




<pK*° 


1 1 K+4.5+I.M 
ix -°-3.7-1.7 


11.11}^ ±0.8 


iu - u -1.5-0.8 


10.7 ± 1.1 


B+ - 


<pK* + 


1U.D_ 49 _ 1 6 


i2.ii^ ± i.i 


a 7+2.1+0.7 
°- ' -1.9-1.0 


9.4 ± 1.6 


B + - 


uK + 


3.21^ ± 0.8 


5.0 ± 1.0 ±0.4 


6.7i^ ±0.6 


5.4 ±0.8 


B + - 


U)7T + 


1.3i^ ± 1.4 


5.4 ± 1.0 ±0.5 


5.9ii;| ± 0.6 


5.3 ±0.9 






10.01^ ± 1.4 


5.3t^ ±0.5 


< 7.6 


5.6 ± 1.4 



Table 3: Branching Fraction results for PV final states including K* , <p or u>. 



itself. The bottom left plots show the distributions of — 21n(£/£ max ) for CLEO II and CLEO 
III combined for the Kir and irn modes with non-zero yields. The remaining right plots are 
two examples from Belle and show the AE distributions in the ities signal region for 7r + 7r~ and 
K + ix~ respectively. The results of the fits used to extract the signal yields are also shown. 

In Table |U the asymmetry results are shown from all three collaborations. All the results 
are compatible with zero and still statistical dominated. 

The additional statistics that will be collected in the next years by BABAR and Belle will 
provide exciting results in the physics of charmless decay. 



mode 




CLEO 


BABAR 


Belle 


WA 


B° -> 


K+TT- 


-0.04 ±0.16 ±0.02 


-0.10 ±0.04 ±0.01 


-0.07 ±0.06 ±0.01 


-0.09 ±0.03 


B+ - 




-0.29 ±0.23 ±0.02 


-0.09 ±0.09 ±0.01 


0.23 ±0.111^4 


0.00 ±0.07 




K%+ 


0.18 ±0.24 ±0.02 


-0.05 ±0.08 ±0.01 


n n7 +o.U9+o.oi 

u - u ' -0.08-0.03 


0.01 ±0.06 


B + - 


K°7T 




0.03 ±0.36 ±0.09 




0.03 ±0.37 


B + - 


7r + 7T~ 




0.30 ±0.25 ±0.04 


0.77 ±0.27 ±0.08 


0.51 ±0.19 


B + - 


7T+7T 




-0.03l£$ ±0.02 


-0.14±0.24i^ 


-0.07 ±0.15 


B + - 


P~7T + 




-0.11+^^0.04 




-0.11 ±0.17 


B + - 


P + 7T~ 




-0.62+^| ±0.06 




-0.62 ±0.29 


B + - 


K + p~ 




0.28 ±0.17 ±0.08 




0.28 ±0.19 



Table 4: Direct CP asymmetry results for Kir, tttc, pir and pK modes 
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